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Abstract 

This paper is devoted to the analysis of the classical Keller-Segel system over M. d , d > 3. We 
describe as much as possible the dynamics of the system characterized by various criteria, both 
in the parabolic-elliptic case and in the fully parabolic case. The main results when dealing with 
the parabolic-elliptic case are: local existence without smallness assumption on the initial density, 
global existence under an improved smallness condition and comparison of blow-up criteria. A 
new concentration phenomenon criteria for the fully parabolic case is also given. The analysis 
is completed by a visualization tool based on the reduction of the parabolic-elliptic system to a 
finite-dimensional dynamical system of gradient flow type, sharing features similar to the infinite- 
dimensional system. 

Key words. Chemotaxis, parabolic systems, global weak solutions, local weak solutions, blow-up, 
energy methods. 
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1 Introduction 

This paper aims to describe the dynamics of the Keller-Segel system in the whole space M. d and high 
dimension d > 3 , with a particular emphasis given to blow-up and related facts. In its simplest and 
nondimensionalized formulation, the system reads as following 

dtn = An — V • (nVc) , 

(KS) 

e dtc = Ac + n — ac, 

and describes, at the macroscopic scale, a population of cells of density n which attract themselves by 
secreting a diffusive chemical signal of concentration c. The nonnegative parameter e is proportional to 
the ratio between the two diffusion coefficients of n and c, appearing in the dimensionalized formulation 
of (|KS|) . It takes into account the different time scales of the two diffusion processes, e = describing 
the chemical concentration evolution in a quasi-stationary approximation. In other words, e is the 
parameter that influences the possible dynamics of c and consequently of n. The chemical degradation 
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rate a is also a nonnegative constant. It is related to the range of action of the signal, as we shall see 
later. 

In the quasi-stationary parabolic-elliptic case (e = 0), the chemical concentration should be under- 
stood as: 

Ed* n , a = , 

(1.1) 

B% * n , a > , 

where hereafter * denotes the convolution with respect to the space variable and Ed and B°[ are 
respectively the Green's function for the Poisson's equation in M. d , and the Bessel kernel: 

Ed{x) =» d W=^ = (d-2)|S*-i| ' (L2) 

r+oo 1 , ,2 

B?w = /„ (L3) 

The system is equipped with non negative initial data no and cq (if e > 0) and with fast decay 
conditions at infinity. Moreover, the initial cell density no is supposed to be an integrable function, so 
that the total initial mass of cells M is conserved along time: 

M = I n (x) dx = n(x,t) dx , 

JR d JM. d 

(see Theorem 11.11 for the complete set of assumptions) . 

The question of blow-up (or concentration) of solutions of the (|KSj) system has been a challenging 
issue in the field of mathematical biology since the seminal papers by Keller and Segel |34t 135] . as 
it describes the propensity of cells to aggregate when they interact through a chemical long-range 
signaling [431 132"! |4"T| I23j . The subtlety of this system has triggered a lot of theoretical works and 
many problems are only partially solved or still open. We refer to |3U| 1431 EZ] for a complete overview 
of various mathematical results and modeling approaches. 

It is well known, since the work by Childress and Percus [16] . that the behaviour of the Keller-Segel 
system strongly depends on the space dimension. In dimension one, blow-up cannot occur for the 
solutions of (jKSp . |40j . In dimension two and for the parabolic-elliptic case (e = 0) we have the so- 
called "critical mass phenomenon": when a = the L 1 — norm of the cell density determines whether 
blow-up occurs or not, i.e. the cell density becomes unbounded in finite time if the mass is above some 
threshold (equals to 8ir), and it does not if the mass is below the same threshold. The same holds 
true when a > up to an additional smallness condition on the second moment of no for blow-up 
only, P3]. This remarkable dichotomy is not yet proved to hold true when e > 0. However, the mass 
still continue to be a key quantity. More precisely, when e > one can prove global existence if the 
mass is below 8ir |13] (using energy methods), but no clear result exists on the other side. In the case 
of a bounded domain, Horstmann proved that blow-up occurs for some initial conditions having very 
negative energy |29|, 128]. To complicate the picture of the parabolic-parabolic case, the authors in [S] 
showed the existence of positive forward self-similar solutions of (|KSp with a = 0, decaying to zero at 
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infinity and having mass larger than 8ir (which is no longer possible in the parabolic-elliptic case), see 
also [12] and the references therein. 

In higher dimensions, the L 1 — norm of no is no longer a critical quantity. One available criterion for 
global existence involves the L d / 2 — norm of the initial cell density |18 t 119 1 f20 | I36|. On the other hand, 
blow-up is shown to occur when e = if the second momentum of no is small enough compared to the 
mass. Nothing is known about the possible blow-up of solutions when e > 0, as in the two-dimensional 
case. 

The main purpose of this paper is to shed a new light on the dynamics of the (JKS]) system in both the 
parabolic-elliptic case and the parabolic-parabolic case, when d > 3. We also propose a visualization 
tool to get the picture of the dynamics in the case e = 0. The tool is based on the reduction of 
(|KSp firstly to a one-space dimensional system and secondly to a finite-dimensional dynamical system 
sharing analogous features with (|KSp . 

More specifically, our main result concerning existence of solutions of (jKSp is the following local 
existence result for the parabolic-elliptic case, without smallness condition on the initial density n$ 
and where the obtained weak solution is not to be intended in the integral sense (mild solution). For 
analogous results concerning mild solutions we refer to |19|. [36] and the references therein. 

Theorem 1.1 (Local existence). Let d > 3, a > and e = 0. Let no he a nonnegative initial density 
in (L l nL a )(M. d ), with a > d/2. Assume in addition that no G L 1 (W d ,ip(x)dx) where ip is a nonnegative 
function such that ip(x) — > +oo uniformly as \x\ — > +oo, e~^ G L 1 (W i ) and |VV>| G L oc (M d ). Then, 
there exists T max > depending only on ||wo||L a (R d ) ana " a wea k nonnegative solution in the distribution 
sense (n, c) of (jKSp . where c is given by (jl.ip . such that 

n€L co ((0,T max );(L 1 nL a )(R d )) and n G L£ c ((0, T max ); U> (R d )) , pG(a,oo). 

Moreover, n G L°°((0, T max ); L 1 (M d , ip{x) dx)), according to the hypothesis on no- In addition, the 
total mass is conserved, 

M= no(x) dx = / n(x,t) dx , a.e. t G (0, T max ) , 

JR d JR d 

nc and n log n belong to L°°((0, T max ); L 1 (E d )), the energy £[n] G L^ c (0,T max ), while the dissipation 
of energy j Rd n|V(logn - c)| 2 G L\ oc (0,T max ). Finally, ifT max < oo then lim \\n(t)\\ L a^ = +oo. 

The above local existence becomes global if in addition we assume that the L d l 2 —norm of the initial 
cell density is small enough. Our contribution in that direction is to obtain a smallness condition for 
l^oll^d Rd improved with respect to the existing one (see Corollary 13.11 and Remark 13. 4|) . Unfortu- 
nately, the results given in Theorem 1 1 . 1 1 and Corollary 13.11 can not be extended to the fully parabolic 
case (e > 0), essentially because in order to do so one has to control the temporal derivative of the 
chemical concentration c (see Remark 13. 5p . 

Concerning the blow-up issue, one of the key ingredient used here is the free energy naturally 
associated to (jKSp . in opposition with the two dimensional case where the free energy is fundamental 
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for proving global existence results only |13| [TU] . The energy functional reads as follows 

£[n,c](t) = I n(x, t) log n(x, t) dx — f n(x,t)c(x,t) dx + - J | Vc(x, t)\ 2 dx + ^ A c 2 (x,t)dx, 

(1.4) 

and it satisfies the dissipation equation 

—£[n,c](t) = -[ n(x,t)\V {log n(x,t) - c{x,t))\ 2 dx - e [ \d t c{x, t)\ 2 dx . (1.5) 

The novelty in this paper, is that we will use also the corrected energy, defined as 

T[n,c](t) = log (J^\x\ 2 n(x,t)dx\ +^£[ n > c ]- ( L6 ) 

Under the quasi-stationary assumption e = 0, the free energy (|1.4p reduces to the difference between 
the entropy and the potential energy: 

£[n](t) = / n(x, t) log n(x, t) dx — — / n(x,t)c(x,t) dx . (1-7) 



2 . 

The blow-up results analyzed here with the help of the above tools can be summarized as follows. 
In the parabolic-elliptic case (e = and a > 0), blow-up occurs in finite time if one of the following 
criteria is fullfilled: 



or 



f |x| 2 n (2;) dx < K?(d, M)M^ , 
I \x\ 2 n Q {x)dx < K 2 {d)M 1+1 d exp (-J-^K 



where the constants Kf(d, M) and K 2 (d) are defined in Propositions 14. ll and 14.21 respectively. After 
having derived both criteria, we show that none of these two criteria contain the other. However 
we have some evidence that the second one (involving the free energy) appears to be better and the 
natural one in some sense. Indeed, this is clearly the case when dealing with the discrete model in 
Section O Moreover, it is possible to extend the second criteria to the parabolic-parabolic case (e > 
and a = for the sake of simplicity). Doing that, we are able to prove the following concentration 
result for the critical norm of the cell density (recovering the parabolic-elliptic case as e \ 0). 

Theorem 1.2 (Parabolic concentration). Let e > 0, a = and d > 3. Assume that the nonnegative 
initial densities (no,co) have finite energy £[no,co] and satisfy 



j \x\ 2 n (x)dx < K 2 (d) M 1+ d exp ( -^rj£[n , c ] ] exp (-e 1 
JRd \ dM J 



(1.8) 



where the constant K 2 (d) is given in (14, 7D and 7 G (0,1). Let (n,c) be a sufficently smooth solution 
of (jKSp generated by (n-o,co) and T max the maximal time of existence (possibly infinite). Then there 
exists a constant C(d) > such that 

sup ||n(*)|| - >^Z^ e -r-i. (1.9) 
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The plan of the paper is as follows. We first list in Section [5] some useful sharp functional inequalities, 
which appear to be key tools for the quantification of criteria involving critical quantities. In Section [3] 
we present the new local in time existence result for the parabolic-elliptic Keller-Segel system without 
any smallness condition on the initial data ( Theorem I l.ip and we review conditions for global existence 
(Corollary 13. In Section 2] we derive the above two blow-up criteria for the parabolic-elliptic case 
(Propositions l4~Tl and l4~2j) and we discuss the complementarity of those criteria (Section I4.2p . In 
Section [5] we extend the criterion involving the free energy to the parabolic-parabolic case into a 
concentration result (Theorem II -2|) . Finally, in Section [61 we shall analyze finite-dimensional energy 
functionals whose gradient flow possesses a structure similar to the parabolic-elliptic (IKS|) system. We 
derive accordingly one criterion ensuring global existence and two criteria for blow-up. The overall 
dynamical picture is drawn in Fig. [TJ 

2 Functional inequalities and preliminaries 

When analyzing the Keller-Segel system (jKSp over M? (whatever is e > 0), the free energy (|1.4p is 
the key tool to prove the global existence of weak free energy solutions (see [9] for the definition). 
Indeed, dual functional inequalities (the logarithmic Hardy-Littlewood-Sobolev inequality and the 
Moser-Trudinger-Onofri inequality) applied to the energy, together with the time decreasing behaviour 
of the energy itself, give the key a priori estimates on the entropy term J R2 ^logn, the potential 
interaction term f^ 2 nc and on the H 1 norm of c, [13] . 

Conversely, when the (jKSp system is analyzed over M. d , d > 3, these inequalities are no longer 
adapted to the energy functional. Therefore, one has to establish a priori estimates on some LP norms 
of the solutions in order to get existence results of more classical weak solutions in the distribution 
sense. 

We shall show in this section that, even in the high dimensional case, when e = sharp functional 
inequalities, corresponding in some sense to those mentioned above for the two dimensional case, are 
the key tools for both the local and the global existence results. More precisely, our purpose is to use 
optimally the sharpness of these inequalities and their "best constant" in order to quantify possible 
thresholds of eventually critical norms of the cell density n. Consequently, the results given here for 
the parabolic-elliptic Keller-Segel system (e = 0) are not only refinement of known results but also 
new ones. This is the case for instance for the local existence of weak solutions without any smallness 
condition on the initial density no and the hypercontractivity property of the solutions. On the other 
hand, the already known smallness condition on no sufficient for the global existence will be improved 
thanks to a more careful estimate and conjectured to be the optimal one (see Remark 13. 4p . 

Concerning the parabolic-parabolic (jKSp system (e > 0), it appears that the estimates valid for the 
parabolic-elliptic case can not be reproduced optimally. Therefore, we are not able to obtain here a 
local existence result without any smallness conditions, nor to improve the known smallness condition 
for the global existence when e > (see Remark 13. 5p . This is one major difference with the two 
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dimensional case, where the energy (|1.4|) can be handle in the same way, whatever is e > 

The mentioned inequalities used here and throughout this paper are firstly the classical Sobolev 
inequality 

<C s (d) II V/IU/Rcn, Cl(d):=- 4 , ■ „,., , (2.1) 



and the following special case of the sharp Hardy-Littlewood-Sobolev inequality [39] 

// f( x )\ x - y\~ X g{y) dxdy < C H Ls(d, A)||/||LP( R d)||p|| LP ( M d) , p = 2d , 0<A<cZ, 
J Jm.dxm.d- 2d -A 

(2.2) 

where the best constant CnLs(d, A) has been obtained by Lieb [38]. More specifically, because of (|1.2p 
and (|1.3p . inequality (|2.2[) will be used in the sequel with A = d — 2. In this case (|2.ip and (|2.2p turn 
to be dual, as proved for instance in |39] Theorem 8.3, with 

so that the following relation can be established 

C 2 s (d) = fi d C HLS (d,d-2) . (2.4) 

Let us observe that when A \ 0, (I2.2p boils down to the logarithmic Hardy-Littlewood-Sobolev 
inequality [15j E] used in the two dimensional case [TOj [13] . 

The other functional inequality we shall use is simply the following Gagliardo-Nirenberg interpolation 
inequality, 

HI 2 (p+d <C GN (p,d)\\Vv\\^ m \\v\\^ , l<p<d. (2.5) 

L p (R d ) { ' LpT" 



We will focus on the particular case p = | in (|2,5p and pay much attention to the best constant 
CGN(d/2, d). Generally, it is not possible to compute the optimal constants in the family of inequalities 
(|2,5p . However and very fortunately, it is possible in the special case p = i under interest here, since 
it corresponds to the Gagliardo-Nirenberg interpolation inequality 

II ,, II 2<t+2 < n 2<r+2\\ v || ad || |,2+(7(2-d) n 2 

for cr = t. Weinstein has characterized the sharp constant in this context [49] . so that 

d 

f d \ ( 2\ 2(3+2) 2_ 

(-,d = C s>d =(l + 5 ) W^), (2-7) 



equation — 4-0 + ipd +1 = 0. 



where tp is a i? 1 (M ci ) n C°°(R ) positive and radial function, solution of minimal L 2 norm of the 



Next come two interpolation lemmas, useful in the sequel for the control of both the entropy 
J Rd n log n and the potential J Rd nc in the free energy. 
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Lemma 2.1 (Entropy lower bound). Let f be any nonnegative L l {U. d ) function such that 
I = J R d \x\ 2 f(x) dx < oo and j Rd f log / < oo . Let M = f Rd f. Then, 

[ flogf + SI> MlogM + ^log(5vr- 1 ), (2.8) 

for any 5 > 0, and 

J /log/ + ^[log/ + 1]> MlogM + lf log (^) . (2.9) 



Proof. Applying the Jensen's inequality with the probability density fi(x) = 5 d ^ 2 ir d l 2 e ^'l 2 , one 
obtain: 

f /log/= I -log(-) / u+ / /log/i 
> MlogM + / /log/., 

and (|2.8|) follows. The optimization of (|2.8|) with respect to 5 > under the fixed constrains / and 
M, yields ([22]), the optimal 6 being 5 = §f. □ 

Remark 2.2 It is worth noticing that the entropy lower bound (|2.8p . and so (|2.9p too, is actually a 
sharp inequality since for any fixed 5 > equality holds true in (|2.8f) if (and only if) / = Mfj>. 

Lemma 2.3 (Potential confinement). Let f be any nonnegative function such that f £ (L 1 C\L2)(M. d ) 
and I = f Rd \x\ 2 f(x) dx < oo. Let M = f Rd f(x)dx. Then, 

2 1 -ImI +1 L 1 -I < [[ f(x)- 1 f(y)dxdy<C HLS (d,d-2)M\\f\\ , (2.10) 



Proof. The right inequality is a direct consequence of the Hardy-Littlewood-Sobolev inequality (|2.2p 
and of standard interpolation (see also [8]). For the left inequality we have simply [1] 



M 2 



< 



f(x)f(y) dxdy 



f(x)\x - y\ 2 f(y)dxdy 



l-2/d 



f(x) 



1 



\x-y 



d-2 



f{y) dxdy 



2/d 



21 M - 2 



xf(x) dx 



l-2/d 



and the proof follows. 



□ 
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We conclude this section by recalling some few basic properties satisfied by the Bessel kernel B% 
defined in (jl.3p and usefull in the sequel. Probably they are classical ones. However, since we do not 
found any references for them, we list and prove the required properties in the lemma below. 

Lemma 2.4 (Properties of the Bessel kernel). The following relations for and V£>^ hold true for 
a > in any space dimension d > 3. 

(i) Expansion of B°[ with respect to E d : 

B2 = E d -aE d *B2 a.e., (2.11) 

(ii) Gradient formula: 

VB2(x) = -^ ZT ^ d g a (\x\), (2.12) 

where g a (\x\) : = is a positive radial function exponentially de- 

creasing from 1 to as \ x | — )• oo. 

(Hi) Corrected Euler's homogeneous function theorem: 

x ■ VB a d (x) = -{d- 2)B2(x) - 2a (B a d * B a d )(x) . (2.13) 

Proof. First of all, let us observe that the convolutions in (|2.1ip and (|2.13p are well defined since both 

d d 

B^ and E d belong to L^ 2 (R d ), (the weak space), and B° € LP(R d ) for 1 < p < (see [22 )■ 
Then, identity (I2.1ip make sense a.e. . Next, identity (|2.1ip is an immediate consequence of the facts 
that (i) E d and B°[ are the Green's functions for —A and (—A + a), respectively, and (ii) B^ * v is 
the unique solution of (—A + a)u = v in V'(M, d ), belonging to L r (K. d ) for some r > 1, (see |39|). 

Formula (|2.12p is a straightforward computation. To prove (|2.13|) . first notice that the Fourier 
transform of B^ is (a + 47r 2 |£| 2 ) , when defining the Fourier transform as /(£) = J Kd e~ 2mx '^f(x) dx, 
(see [39]). Then, denoting / v the inverse Fourier transform, we have in the Fourier variable £, 

(x ■ V x B2r(0 = (x ■ V x ((a + 47r 2 |e| 2 )- 1 ) v )^(0 = 27T, (x ■ (£(a + 4^|£| 2 r 1 ) v H£) 



a + 47r 2 \^J~ a + 4vr 2 |e| 2 ( a + 4^ 2 |£| 2 ) 2 
d-2 2a 



a + 4vr 2 |e| 2 (a + 4 7 r 2 | 1 e| 2 ) 2 
and identity (|2.13p is proved. 

□ 

Remark 2.5 Actually, in order to prove Lemma [2.4[ it is possible to develop an alternative argument 
based on the homogeneity property of B^, namely: B^iXx) = \ 2 ~ d B2^ ' (%)■ From this point of view, 
formula (|2.13|) can be seen as a corrected Euler's homogeneous function theorem. Obviously (|2.13p 
reduces exactly to the Euler's homogeneous function theorem for E d as a = 0, E d being homogeneous 
of degree (2 — d). 
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3 Local and global existence for the parabolic-elliptic system (e = 0) 



This section is entirely devoted to the proof of Theorem 11.11 The proof is given in several steps, 
beginning with a cascade of a priori estimates. The functional inequalities given in the previous 
section and the mass conservation will play a fundamental role. The rate of decay at infinity of 
the initial density no is also essential. Let us observe that one can weaken the usual hypothesis 
n £ L l (R d , \x\ 2 dx), as we have done here, and obtain existence result also with initial density having 
infinite second moment but logarithmic decay (see also Subsection 13.51 and Remark 13. 3p . This is 
absolutely not the case when dealing with blow-up results, always obtained under the hypothesis of 
finite initial second moment. 

As a byproduct of the a priori estimates, we obtain the following global existence result under an im- 
proved smallness condition on the initial density. Since the corollary is a straightforward consequence 
of Theorem 11.11 the proof will be omitted. 



Corollary 3.1 (Global existence). Under the same hypothesis on the initial density uq as in Theo- 
rem li.il if in addition 

8 2(1 + |) (d 



l«oll,*™< ~, C GN d '[^d) , (3.1) 



the above weak solution exists globally in time. Moreover, the Lz-norm of n is time decreasing. 



Remark 3.2 In high dimension (d > 6), Corollary 13.11 can be proved under the sharp hypothesis on 
the initial density no £ {L 1 r\L^){M. d ), as it follows by an attentive analysis of the proof of Theorem 11.11 
see subsection! 



3.1 Local in time estimate of ||n||x,p( K d) 

Since —Ac = n — ac , we have for any nonnegative solutions n 

^- [ n P = -^l [ |Vni| 2 + (p-l) f n p+1 -a(p-l) [ n p c . (3.2 
Next, using standard interpolation and the Sobolev inequality ()2. 1 j) . we get 

for any max{| — 1 ; 1} < p < oo. Inequality (|3.3p gives us for p > | 

where r = ^j, r' = and 5 > 0. Plugging estimate (I3.4p into (13.2f) . one obtains 



(3.3) 



nP < (p-1) { dCUd) 



P 



26 r 



|Vn5||^ 2(Rd) + (p-l)_( jy) (3.5) 



i+Cp-I)- 1 
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Then, for any fixed p > i, it is enough to choose 5 = 5(p) > such that ^fr*^ — 4 < in order to 
have from (|3.5p a local in time boundeness of the IP norm of n, p > |, whenever no E L p (R d ), i.e. 

n p (t) < h p (t) [ n p , (3.6) 



for t S [0, T p ), where 

(p-I)- 1 



(f-p) 



(I-p)- 



(3.7) 

This is obviously only a local in time estimate since h p (t) — > +oo as t /*■ T p . 
3.2 Global in time estimate of ||ri||ii>nRd> under smallness condition 



In order to obtain a global in time estimate of ||n||£ P ( R d), one has to be more careful in the estimate 
of the L p+1 norm of n than in (|3.3p and use directly the sharp Gagliardo-Nirenberg inequality (|2.5p . 
Indeed, we have 

L n,+1 = L (n ' )lir+1> - c ^ 1,(M) l|VnSfe « M w<> • (3 - 8) 

Plugging (|3.8|> into fl3.2|> we get for any 1 < p < d 



^ /" < (P- l)[cl^ +1) (p,d)|h(t)|| i# ^ - 1] HVrxf ||| 2fRd) . (3.9) 



dt 



P 



As a consequence, choosingp = | into (|3.9p . whenever the initial density satisfies (|3.ip . the ||w(f)|| d 
norm decreases for all times t > 0. 

3.3 Local in time hyper-contractivity property 

Let no £ (L 1 D L a )(R d ), with a > | arbitrarily closed to i. From subsection 13.11 there exists a finite 
time T a > 0, defined in f[3"7f|) . and a local solution n G L°°((0,T a ); (L 1 n L a )(M d )). We are going to 
prove that, for any p S (a, +oo), there exists a constant C not depending on ||wo||LP(M d ) such that 



f n p (t) < C(l + t 1-p ) , a.e. t £ (0,T a ) . (3.10) 

Let d € (0,T a ) be arbitrarily small. Owing to the local in time boundedness of the L a -norm of n, 
there exists a modulus of " ^-equintegrability" u)(K\ T a — 5) such that for K > 1 it holds 

sup \\(n- K) + {t)\\\ <uj{K;T a -5) and lim u(K; T a - S) = . (3.11) 

0<t<(T a -6) L"3 (R d ) A'->+oo 

Indeed, for £ S [0, T — 5], using the increasing behaviour of h a (t) defined in (|3.7p . we have 
\\(n-K) + (t)P <* [ n a (t)< ha{Ta ~ 5) [ nt =: lo(K; T a — 5) . 
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Next, let us estimate the time evolution of \\(n — K) + \\ LP ^ Rd ^, from the differential inequality 
| J Rd {n - Kf + < -4^il f Rd |V(n - K)\\ 2 + (p - 1) J Rd (n - 

(3.12) 

+K(2p - 1) / Rd (n - K)l + p/<T 2 J Rd (n - Kf' 1 . 

For any fixed p > max{2; a}, interpolating the LP and L p_1 norms of (n — between L 1 and L p+1 , 
one get easily from ()3. 12[) 

4- [ (n-K) p + < -4 (P ~ ^ [ \V{n-K)l\ 2 + 4(p-l) [ {n - K) p + l + AK P M . (3.13) 

dt JR d P JR d JR d 

It remains to take advantage of the negative term in the right hand side of (|3.13p . From standard 
interpolation and the Sobolev inequality (|2.ip we have for any p > max{f - 1; 1} 

v p+1 < \H P vd \\v\\ d, = \\v%\\ 2 2d \\v\\ d, A < <?|(cf) ||V^|| 2 \\v\\ ir „ N . (3.14) 



Let us observe that here it is more convinient to use (|3.14p than (|3.8p . in order to have a larger range 
of p index for which the inequality holds true. 
Using (|3.14p and (|3.1ip . we obtain for t £ [0,T a — 5] the estimate 

/ (n - K)l +1 < C 2 (d)\\V(n - K)\\\ 2 d) u>$(K;T a - 5) , 

JR d v ; 

so that (I3.13P becomes 

f (n-K) p + <-4(p-l)\p- 1 Cs 2 (d)iO-^(K;T a -5)-l\[ (n — K) 1 ^ 1 + \K V M . (3.15) 



eft 



Taking K sufficiently large, the quantity rj := p 1 C s 2 (d)uj " d {K\T a — 5) — 1 



is positive and using 



once again the interpolation of the LP norm of (n — K) + between L 1 and L p+ , we obtain 

^- I {n-K)\ < -4(p- 1)t/M _ J^ ( /" ( n _ jqp V _1 + 4 K p M . (3.16) 

jR d \JR d J 

We are finally able to prove that there exists a positive finite constant C = C(T a ,K,M,p), not 
depending on / K d(^o — ^0+> such that for any p > max{2; a} and for t G (0, T a — J] 



simply by comparison of positive solutions of (|3.16j) with positive solutions of the differential equation 

i p 

v!(t) +4(p— 1) ?7 M p- 1 up- 1 (t) = AK P M . Consequently and as usual, the hypercontractivity estimate 
(|3.10p holds true for any p > max{2; a}. For p G (a, max{2; a}] it follows by interpolation. 
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3.4 A priori estimates for the chemical density 

Let n G (L 1 n L a )(]R ci ). Then, from the previous estimates, any corresponding nonnegative solutions 
n belongs to any L p (R d ), p G (l,oo), for a.e. t G (0, T a ). Consequently, taking p G (1, ^), c in (ll.ip is 
well defined for any a > 0, it belongs to L q (R d ), with q = , and 



L9 



< C(d,a) \\n(t)\\ LP(Rd) , a.e. t G (0,T a ) , 



by the weak Young inequality [39], (see also Lemma l2.4p . In the same way, Vc is well defined in the 
distributional sense as 

f VE d *n, a = 0, 
Vc= { 

[ VB% * n , a > . 

and |Vc| G L r (R d ), with r = for p G (1, d): 

\\Vc(t)\\ L r m <C(d,a)\\n(t)\\ LP{Rd) , a.e. te (0,T o ). (3.17) 

3.5 Entropy, potential, energy and energy dissipation estimates 

d 

To conclude the a priori estimates, we will show how a local in time control on the La norm of n 
together with a control on the decay of tlq as |x| — > oo, give local in time estimates on the entropy, 
the potential, the energy and the dissipation of the energy. 

Recall that under the assumptions of Theorem 11.11 we have that the initial density satisfies no G 
{L l C\Li )(M d )r\L 1 (M d , Tp(x)dx), where ip is nonnegative and satisfies e - ^ G L 1 (R d ) and |V^| G L°°(R d ). 
Let n be a nonneg ative solution such that n(t) G (L 1 n L2)(R d ) for a.e. t G (0,T), T > 0. Then, the 
local in time bound on the potential follows directly from the definition of E d and (|2.1ip (for the case 
a > 0) and from the potential confinement Lemma 12.31 that gives 

0< / n(t)c(tt)<fx d C HLS (d,d-2)M\\n(t)\\ a (3.18) 

J R d L?(K a ) 

Next, the positive and negative contributions of the inital entropy J Rd no log no are bounded. Indeed, 
for the positive contribution it simply holds true, for any d > 3, that 

0< / (no log nn)+ < linn II d , . (3.19) 

Moreover, setting v = no l{ no <i} and m = J* Rd the Jensen inequality gives us 

(v log v - vip) = / n log (—r) > mlogm- mlogdle^H^i^dO, 

which implies 

0< / (n logn )- < mlog(||e - ^|| i i/ R d))-mlogm- / n l {no <i)^ < C(m, ||e" l/ '|| i i (IR[ i ) )+ / n ^ . 

JRd JRd JRd 

(3.20) 
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The previous computations (|3.18j) . 19|) and (|3.2U|) give us that the initial energy (|1.7|) is bounded as 
follows 

-C(m, He^Hii^d)) - / n ijj - C {M , d)\\n \\ ^ < S[n ] < \\n \\ L d 



Now, let us extend the above estimates locally in time. Estimate f|3. 19j) is obviously true for n(t). 
Concerning the weighted L 1 -norm of n we have 



— f nil) = — f nVtfj ■ V(log n — c) < — f \Vifj\ 2 n + 5 f n|V(logn — c 
with arbitrary 5 > to be choosen later. Hence, it follows that 

/ n(t)i/>< [ no4 + h\Vi/>\\ 2 Loo(R *) Mt + S [ I n|V(logn-c)| 2 . (3.21) 
It remains to estimate the energy dissipation from (jl.5p . Using (|3. 18j> . (|3.20p and (jl.7p . we obtain 
Jo Sm d re ( s )l v (log n(s) - c(s))\ 2 dx ds = £[n ] - £[n](t) 

- N^o 11^ d (IRrf) + / H<I (n(t) log n(t))_ + \ J Rd n(t)c(t) ( 3-2 2) 
< ho\\d +C{m,\\e-^ Lim ) + f m MW + C{M,d)\\n{t)\\ * 



Plugging now (|3.2ip with any 5 G (0, 1) into (|3.22p . we finally get 

(1 - 5) J* j n|V(logn - c)| 2 < C(n , t/,, M)(l + t) + C(M, d)||n(t)|| £ $ ^ . (3.23) 
Consequently, from (|3.2ip . the same estimate holds for the weighted L 1 -norm of n , i.e. 

n(t)i> < C(no, iP, M)(l + 1) + C(M, d)|Kt)||^ f ^ (3.24) 
and this gives us the control on the energy 

- C(n ,V,M)(l +t) - C(M,d)||n(t)|| , < £[n](t) < ||n(i)||^ 



and exactly the same control on the entropy. 

Remark 3.3 (Logarithmic decay of the density) The same estimates of this subsection hold true if 
we assume that no G L 1 (M' i , |x| 2 <ix) but with slightly modified computations (see [lQj[13]). Using here 
the hypothesis no G L 1 (M. d , ij)(x)dx) we want to underline that one can obtain local and consequently 
global existence results under weaker condition on the decay at infinity of the initial density. For 
instance, it is sufficient to consider logarithmic decay at infinity for no, i.e. ip(x) = dlog(l + |x| 2 ). 



13 



3.6 Proof of Theorem 11.11 

To prove Theorem 11.11 we need to regularize system (|KS|) with e = and c given by (jl.ip and to 
show that the sequence of regularized solutions satisfies firstly the previous estimates and secondly is 
relatively compact in an appropriate topological space. This is now a quite standard procedure, used 
in the context of Keller-Segel system for example in [131 EE EO]- However, we give here a quite 
rapid sketch of the proof for the sake of completeness. 
As a regularized problem we consider 

d t n a = An ff - V • (n a Vc a ) , (3.26) 

where c a is given by 

[ {E d *n°(t)*p°)(x), a = 0, 
c a (x,t) = < (3.27) 
[ {B% *n a (t)* p a )(x) , a>0, 

and p a is some sequence of smooth positive mollifiers with Hp^H^i^d) = 1. The regularized initial 
condition 

nl = n * p a (3.28) 

is also considered. Problem (|3.26p - (|3.27p - (|3.28p has a nonnegative smooth solution as follows by the 
Schauder's fixed-point theorem. Moreover, the solution (n CT , c u ) satisfies all the a priori estimates given 
in the previous sections. Indeed, one has essentially to check that the master equation (|3.2p holds true 
at least as an inequality, so that the fundamental estimates (|3.3p and (|3.8p on the L p+1 norm of n a 
can be applied. This is the case since 

-/ (n a ) p Ac a < [ (n a ) p (n CT * p a ) < [ (n a ) p+1 . 

Concerning the compactness of the sequence {n a }, we are intended to use the Aubin compactness 
lemma in [2]. Therefore, we chose B = L 2 (R d ), X = H l {R d ) n L 2 (R d , ip(x)^dx) compactly imbedded 
in B and Y = H^ 1 (R d ) so that B C Y. Using the previous a priori estimates, we prove below that 
{n ff } is bounded in L 2 ((5,T a — 5); X) uniformly in a, where 5 £ (0,T o ) is arbitrarily small and T a is 
defined in (|3.7p . The boundedness of {9tn CT } in L 2 ((5,T a — 5);Y) uniformly in a follows then from 
the same computation. Consequently, the Aubin's lemma gives that {n a } is relatively compact in 
L 2 ((5, T a — 5); B) and the proof is complete. Let us observe that for high dimensions (d > 6) we obtain 
in fact that {n a } is relatively compact in L 2 ((0, T a — 5);B). 

For the sake of simplicity, we omit the index a in the sequel. First of all, let us observe that 
n G L°°((5, T a -S); L 2 (R d )) if d = 3 and n G L°°((0, T a -S); L 2 (R d )) if d > 4, as it follows by (l3~T0|) 
and by (|3.6|) respectively, being a > |. 

Next, from (l3~T7l) with p = f, we have that |Vc| G L°°((0,r a - 5) ; L d (R d )) , for any d > 3. 
Consequently, since by the Holder's inequality, 

||(nVc)(t)|| L 2 (R d) < \\n(t)\\ l _2^ \\Vc(t)\\ Ld{Rd) , 
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using again (EHUD and QMS) we obtain that ra|Vc| G L°°((5,T a - 5);L 2 (R d )) for d = 3,4,5 and 
n|Vc| G L°°((0, T a - 5) ; L 2 (M d )) for d > 6. 

Finally, multiplying the equation on n in (|KSp against n, integrating over M rf and then over (5, T a — 5), 
one easily obtain 

rTa—6 rTa—6 

/ \\Vn(s)\\ 2 L2{Rd) ds <\\n{5)\\ 2 L2{Rd) + \\(?iVc)(s)\\ 2 L2{Rd) ds , 

i.e. |Vn| G L 2 ((5,T a -S);L 2 (R d )) for d = 3,4,5 and |Vn| G L 2 ((0,T a - <5);L 2 (M d )) for d > 6. 

It remains to estimate the L 2 (M. d , if)(x) 3 dx) norm of n. This is an immediate consequence of the 
computation 

^™ 2 (t) < ( / ^n(t)V ( / n 3 (t)V , 

and of estimates (pM) and again (15TUD and ([33]) . Therefore, n G L°°((J,T a - 5) ; L 2 (R d ,ip(x)^dx)) 
for d = 3,4,5 and n G L°°((0,T a — S) ; L 2 (M d , -0(x)^dx)) for d > 6. 

□ 

Remark 3.4 (Improved smallness condition) In [20], in order to obtain the global existence result, 
the authors have used inequality (|3.14p instead of (j3.8[) and the corresponding inequality 



d_ 

dt 



n p (t) < (p-l)[cUd)\\n(t)\\ Lim - J]||Vni||i 2(Rd) . (3.29) 



instead of (|3.9p . with the corresponding smallness condition 



|noll ^w < dcM • (3 - 30) 



2(1+ — ) 

Dince we can prove that C G ]v d! (§,d) < C|(d), condition is a weaker condition than (|3,30p . 

The last together with the facts that 

(i) inequality (13.8|) is sharp (equality holds true for the minimizers of (12. 6h ) ; 

(ii) under the smallness condition (|3.30p : all the 1 1 ^C*) 1 1 norms are time decreasing, for p G 
(1, min{p*, d}] and where p* := 4||no|| _ d Cg 2 (d) is greater than |; the hyper-contractivity 

property (|3.10p holds true globally in time for p > min{p* , d} ; the entropy L d n log n is time 
decreasing in dimension d = 3; 



induced to conjecture that the smallness condition (|3.ip is the optimal one. 

To prove the claimed inequality between optimal constants, it is sufficient to write the Gagliardo- 
Nirenberg inequality (|2.6p for a = | and the function v for which the equality holds true, i.e. 

2(1+|), ||2 ii A _ || ,|2(1+|) 
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Observing that by standard interpolation and the Sobolev inequality (|2.ip we have 



2( lH - — ■ ~) 

and using 

flgHD , the inequality C G ( N d> (f , d) < C|(d) is proved. Since equality cannot hold true in 
(|3.14p while it can in (|3.8p . the above claimed strict inequality is obtained. 

The time decreasing behaviour of the family of L p -norms of n, p £ (l,min{p*,d}], whenever the 
smallness condition (|3.30|) holds true, is a direct consequence of inequality (|3.29|) . where the constant 
Cg(d) doesn't depend on p. In order to obtain the same result under the improved smallness condition 

(|3.ip . one should know the behaviour of the constant Cq^ + (p, d) in (|3.9p with respect to p. 

The hyper-contractivity property (|3.10p holds true globally in time for p > min{p* , d} since we can 
choose the modulus of "^-equintegrability" in (13. lip time independent. 

Finally, let us show that (j3.30f) implies the time decreasing behaviour of the entropy J Kd n log n, at 
least in dimension d = 3. Indeed, 

I nlogn = -i[ |Vv^| 2 + / n 2 . (3.31) 
dt JR d Jm d JR d 

Then, using the Gagliardo-Nirenberg inequality (|2.5p with p = 1 and v = y/n to estimate J Kd n 2 in 
(|3.3ip . we have 

j t f^n\ogn< ||V^||i 2(Rd) [c^(l,d)||n|| L4(R(1) -4] . (3.32) 

The best constant CQ^(l,d) is unknown (at our best knowledge) except for d = 3 (see [21]). In that 
case Cgjv(1,3) = (^t 7 ) 1 ^ 6 , while C|(3) = ^(^p) 2 ^ 3 from (12. lj) and the issue is proved since the right 
hand side of ()3.32j) turns to be negative under the smallness condition (|3.30|) . 

Remark 3.5 In the parabolic-parabolic case (e > 0) we can not have the local existence result 
without smallness condition on the initial densities, obtained for the parabolic-elliptic system. This is 
due essentially to the fact that in the master equation (|3.2p one has to add the term — e (p— 1) J" Rd n p dtc 
stronger with respect to J* Rd n p+1 , since the only control that we have for the temporal derivative of 
c is an L 2 control from (|1.5p . For local and global existence results of (IKSp with e > 0, we refer 
to [I9l [37] and the references therein. 



4 Blow-up for the parabolic-elliptic system (e = 0) 

In order to characterize blowing up solutions of the parabolic-elliptic (|KSp system, the general idea is 
to follow the evolution of the second moment of n, i.e. I(t) := J Rd \x\ 2 n(x,t) dx, and to prove that I 
satisfies a differential inequality of type 

U(t) < f {!{€)) , (4.1) 
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where / is a continuous nondecreasing function such that /(0) < 0. Indeed, whenever we can exhibit 
(|4,ip . defining /* := inf{/ > 1 /(/) = 0} G (0, +oo], for any sufficiently smooth solution n of 
(|KSp with finite initial second moment satisfying 1(0) < I*, there exists a time T* < oo such that 

lim I(t) = 0. The vanishing behaviour of the second moment implies the blow-up of some norm 

t/T* 

critical for the existence of the solution and henceforth of the solution itself. Obviously the blowing 
time is in general smaller than T* and at our best knowledge it is an open problem to characterize 
the blowing time in term of the exploding critical norm instead of 1(0). 

The above technique as been applied firstly by Biler in [1] for a model of gravitational interaction of 
particles on a star-shaped domain of M. d , similar to the (|KSp system with e = a = 0. Successively, it 
has been used by several authors in the context of the Keller-Segel system (see for instance [4*01 |2C"1 [TO] ) . 
The methodology is also reminiscent of the blow-up criteria for the nonlinear Schrodinger equation 
initiated by Glassey [21], (see also [M]), and has been successively applied to kinetic gravitational 
models [25] and kinetic chemotaxis models 

Notice that there exists an alternative (and non constructive) method to obtain the existence of 
blowing up solutions, based on energy features (unbounded from below) and the analysis of the 
possible stationary states. We refer to "2H1 EI] for more details. 

Concerning the parabolic-elliptic (|KSp system in two space dimensions, some candidate function / 
in (|4.ip can be explicitly computed, whatever a > is, and the condition /(0) < reads as M > 87T 
when a = 0, (see [101113]). This threshold for the mass M is sharp since global existence can be proved 
for M < 8tt and suitable conditions on the initial density no, ([El 03]). When blow-up occurs, the 
solution ceases to exist classically and possible extensions after the blow-up time have been proposed 
in [46 1 147 ] 122] . depending upon the system regularization. Finally, in the critical case M = 8n, it has 
been proved in [9] that weak free energy solutions still exist globally in time and that they blow up in 
infinite time. 

In the case d > 3 the derivation of (|4.ip is a little more complicated since the potential appears in 
the evolution equation for the moment /. More specifically, from 



d_ 

di 



I(t) = 2dM + 2 n(x,t)x-Vc(x,t)dx, (4.2) 



we have (after symmetrization of the integral term) for a = 



d 1 f f 1 

- I(t) = 2dM - — p- / / n(x, t)- r^niy, t) dxdy (4.3) 

at p i j jR d xM. d \ x ~ y 



while for a > 



5'M = *w - pii, // RW »(«, ^l^nfe, , (4.4) 

where g a is defined in Lemma 12.41 Therefore, different blow-up criteria can be obtained according to 
how the right hand side of (]4.3h or (14. 4p is estimated with respect to /, £[n] and M. Let us observe 
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that (|4.4|) becomes (|4.3p when a = since <7o(M) = 1- The continuity with respect to the dimension d 
is also satisfied in the sense that from (|4.3p and (|4.4p we obtain the evolution equation of / when d = 2. 

In the sequel we shall derive two inequalities of type (|4.ip together with the consequent blow-up 
criteria and we shall discuss their complementarity. One of the criteria is yet contained in |20l and 
both are contained in [J] for the case of star-shaped domains. However, only the case of the parabolic- 
elliptic (IKSp system in absence of chemical degradation (e = a = 0) was considered in the cited papers. 
Here we improve the criterion involving the free energy (not required to be negative) and we extend 
both to the case a > 0. 

4.1 Derivation of two blow-up criteria 

Proposition 4.1 (First blow-up criterion). Let e = 0, a > 0, d > 3 and a > |. Assume that the 
nonnegative initial density no G (L 1 n L a )(M, d ) has finite second momentum satisfying 

\x\ 2 n (x) dx < Kf{d, M)M^ , (4.5) 

where K"(d, M) is defined in (j4.9|) for a > and 

K°(d,M) = K^d) := 2~^(d|S d - 1 |)~^ . (4.6) 

Then, the solution to the parabolic- elliptic (|KSp system constructed in Theorem \l.l\ blows up in finite 
time, that is the maximal time of existence T max is finite and lim ||?i(i)llL a (R d ) = +oo. Moreover, 

t /Tmax 

the blow-up condition J^.5[ ) is incompatible with the smallness condition \3.30\i for global existence. 

Proposition 4.2 (Second blow-up criterion). Let e = 0, a > 0, d > 3 and a > i. Assume that the 
nonnegative initial density uq E (L 1 n L a )(W i ) has finite second momentum satisfying 

[ \x\ 2 n (x)dx <K 2 (d)M 1+ i exp(--^-£[n ]] , K 2 (d) := . (4.7) 

J R d \ dM J 2tt 

Then, the solution to the parabolic- elliptic (|KSp system constructed in Theorem \l.l\ blows up infinite 
time, that is the maximal time of existence T max is finite and lim ||?i-(i)llL a (R d ) = +°°- Moreover, 



the blow-up condition {J^. 7| ) are incompatible with the smallness condition (I3.30P for global existence. 

Proof of Proposition [7^1} Assume first for simplicity that a = 0. Using the left inequality in (12.101) 
into (|4.3p . we have: 

-lit) < 2dM - | S rf-i|-i 2 i-i M -2 +l L l --2(t) , 
dt 

i.e. (|4.ip holds true with the increasing function /(A) = 2dM — |S rf ~ 1 |~ 1 2 x ~ 2 M ? +1 A 1_ 2 satisfying 

d 

f(\) — > —00 as A \ and /(/*) = with /* := K\{d) M d ~ 2 . Hence, we obtain an obstruction to 
global existence when 7(0) < I*, which is the condition (|4.5p . From estimates (|2.9p and ()2.10p . we 
deduce the existence of T\ > 0, T 2 > and T max > such that 

lim / n log n = +00 , lim / nc = +00 and lim ||n(t)|| ra (Md\ = +00 
tSTiJnd t^T 2 J R d tST max " " ( ' 
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with T max < min{T 1 ,r 2 }. 

When q > 0, in order to use the potential confinement Lemma as before into (|4.4p . because of the 
decreasing behaviour of g a , one has to proceed as following 



I s I J J\x-y\<R \ x ~ V\ 

< 2dM - 9 ,^) R }, 1 1 n(x,t)- l —j Z ^n(y,t)dxdy 

l§ I JJ\x~y\>R \ x - Vl 

Then, using the left inequality in (|2.10p . we obtain for any R > 

-I(t) < 2dM - 9 -^- 2 1 ~-2M-2 +1 I 1 -i(t) + , ? a ( R \ 9 M 2 . (4.8) 

The differential inequality (|4.8|) is again of type (|4.ip . with in the r.h.s. an increasing function of /, 
converging to — oo as / \ and to a positive constant as / goes to +oo. To obtain the best blow-up 
condition, one has to maximize with respect to R the value I*(R) for which the r.h.s. of (|4.8p vanishes. 

d 

Such a computation give the definition I* := K®(d, M)M d ~ 2 , where 



Xf(d,M) :=i 



sup 

i?>0 



d-2 



2d|S d - 1 | J R d - 2 + 5a ( J R)M 



2 
d-2 



(4.9) 



Owing to the properties of g a for a > 0, the function to maximise into formula (14. 9p is a positive 
continuous function, vanishing for R = and decaying to as R — > oo. Therefore, the supremum in 
(|4.9p is achieved and iff is well defined. Consequently, the blow-up criterion (|4.5p has been obtained. 
Moreover, Kf(d, M) = K\(d) because go(M) = 1 and the supremum in (I4.9P is achieved for R — > +oo. 

Finally, let us prove the claimed incompatibility. It is sufficient to consider the case a = since iff 
is strictly decreasing with respect to a so that Kf(d,M) < K\{d). From the potential confinement 
Lemma 12.31 and from (14. 5ft . we obtain 

||n<>y > 2 1 -! Mi C- H \ s {d, d-2) /^(O) > 2d\§ d - 1 \Cj I 1 LS (d, d-2). (4.10) 



Using the relation (|2,4p into (|4.10p . the reverse of condition (|3.30p follows. 

□ 



Proof of Proposition \4-%\ Let us again consider firstly the case a = 0. Then, we shall make use of 
the definition of the free energy, its time decreasing behaviour and of the entropy lower bound (|2.9p 
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instead of the potential confinement inequality into (|4.3p . to get 

4-I(t) = 2dM - (d - 2) f n(x,t)c(x,t)dx 
dt J R d 

= 2dM + 2(d-2) (s[n](t)-J n(x,t) log n(x,t) dx) (4.11) 
< d(d - 2)M log I(t) + 2{d - 2)S[n ] + B(d, M) , 
where the constant B(d, M) is defined as 

B(d, M) := d 2 M - 2(d - 2) 



MlogM + — log f — 



(4.12) 



Therefore, (|4.1[) yields with the nondecreasing function /(A) = d(d-2)M log A+2(d-2)£[n ]+-B(d, M), 
satisfying /(A) -)• -oo as A \ and /(/*) = with I* := exp (—^£[n ] - d(d-2)li ) • A § am > for 
1(0) < I*, that is condition (|4.7p . we obtain an obstruction to global existence as before. 

When a > 0, we proceed as for a = 0, but using formula (|2.13p for the Bessel kernel into the 
evolution equation of / write as (|4.2p . Therefore, after the symmetrization of the integral term in 
T2l) . we have 

d 



d -I(t) = 2dM+ll n(x,t)(x-y)-VB2(x-y)n(y,t)dxdy 



2dM - (d - 2) / / n(x, t) B%(x - y) n(y, t) dxdy 
-2a [[ n(x,t)(B% * B%)(x - y)n(y,t) dxdy 

J JR d xR d 

2dM-(d-2) n(x,t)c(x,t)dx -2a (B% * nf(x,t)dx 
Jw. d Jm. d 

2dM + 2(d- 2) ( S[n](t) - \ n(x, t) log n(x, t) dx ) - 2a I c 2 (x,t)dx. 



Neglecting the last negative contribution we are reduced to the previous estimate (|4.1ip and we can 
conclude as above. 

Finally, let us prove the claimed incompatibility. Applying the logarithmic function to both sides 
of (|4,7p and using the definition of the energy, we arrive easily to the inequality: 

dM , , , 4 f , dM , fdM\ If . . . . , dM 

no(x)co(x) dx- 



(logl(0) + l) + J^n {x) log n {x)dx-M log M-^- log (^) < \ 



2 /rod 2 V 2tt ) 2 /rod d — 2 

(4.13) 

The left hand side of (|4.13p is nonnegative owing to Lemma 12.11 Then, from ()4.13|) and again the 
potential confinement Lemma 12.31 we have 

[ <l n (x)c (x)dx < fi d C H Ls(d,d-2)M\\n \\, (4.14) 

JR d L ( K > 



d-2 

and the reverse of condition (|3.30p follows. 



□ 
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Remark 4.3 We are actually not able to prove that both criteria, (|4.5p and (|4.7|) . are also incompatible 
with the improved smallness condition (|3.ip for the global existence. This is due to the lack of 
knowledge about the constant CcN(d/2,d). 



Remark 4.4 (i) It is easy to check the invariance of criterion (|4.5p for a = and of criterion (j4.7|) 
for any a > 0, under the scaling tiq(x) — > tiq(x) = \~ 2 no(\~ 1 x) , preserving the La— norm. This 
invariance doesn't hold true for criterion (|4.5|) when a > because of the dependency on M of the 
constant Kf. However, since Kf(d,M) is strictly decreasing with a, so that Kf{d, M) < K\{d) for 
any a > 0, and because of the claimed incompatibility, it is not possible to construct from a density 
no satisfying the smallness condition (|3.30p a density n^x) = \~ 2 no(X~ 1 x) satisfying criterion ()4.5j) 
for some a > 0. (ii) The decreasing behaviour of Kf with respect to a is naturally due to the fact 
that the degradation term —ac in the (IKSp system prevents blow-up of the density n. This term has 
been completely neglected in the deduction of the second blow-up criterion (|4.7p involving the energy. 
It would be interesting, eventhough hard to do, not to neglect it. 

Remark 4.5 Looking attentively at the proofs of Propositions 14. ll and I4T21 it appears evident that 
the second blow-up criterion is sharper than the first one. Indeed, the first blow-up criterion has 
been obtained with the help of the potential confinement inequality (|2.10p , while the second blow-up 
criterion has been obtained applying the entropy lower bound (|2.9p and that the latter is a sharp 
inequality (see Remark I2.2p . 

Remark 4.6 (Corrected energy) Coming back to the line (|4.11|) and using the fact that j- t £[n(t)} < 
we deduce the following differential inequality for the corrected energy defined in (jl.6p 



In Section[5]we will generalize inequality (|4.15p as well as the blow-up condition (|4.16|) to the parabolic- 
parabolic Keller-Segel system (|KSp in order to obtain a concentration result for the La-norm of n. 

4.2 Complementarity of the blow-up criteria 

The goal of this subsection is to construct examples of initial data no satisfying either 



F[n](t) < d{d-2)MF[n]{t)+B{d,M) 



(4.15) 



Moreover, the blow-up criterion involving the initial free energy (j4.7j) reads equivalently as: 



d{d - 2)MF[n ] + B(d, M) < . 



(4.16) 



K 2 (d)M 




x\ 2 n (x) dx < Ki(d) M^-a , 



(4.17) 



or 




JR d 



x\ 2 n (x)dx < K 2 {d)M l+ d 




(4.18) 
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Here a = 0. The case a > is not considered for sake of simplicity. 

We start with the latter (|4.18p which is more natural, since it is a direct consequence of the fact 
that the energy is unbounded from below and of the mixed homogeneities of the potential and entropy 
terms composing the energy with respect to the mass-preserving dilation f(x) — > / A (x) = \~ d f(\~ l x). 
For this purpose, we introduce the following family of densities indexed by A > 



<{x) := - 



+ A d ip 



where a ^ is some point to be chosen later and (p is a nonnegative function in (L 1 n L? 
that L d ip{z) dz = M and Suppi/? C B(0, 1). Then, the densities Uq belong to (L 1 n L 

d 

mass equal to M, La -norm increasing as A \ and the second moment given by 



(4.19) 

such 
, have 



/ 



\x\ 2 nn(x) dx 



/ \a + \z\ 2 ip(z) dz + / \—a + 



\z\ 2 Lp(z)dz = M\a\ 2 + A 2 / \z\ 2 (p(z)dz. 

(4.20) 

When evaluating the free energy, the cross-interaction between the two densities located around a 
and —a is zero in the entropy term, if A is small enough so that the supports B(a,\) and B(— a, A) 
are disjoints. Hence, we have 



ip(z) log ip(z) dz — dM log A — M log 2 



ip(z) 



1 



z 



/Id— 2 



4 



2a + A(z 



Id— 2 



y(z') dzdz' 



<f(z') dzdz' 



(4.21) 



which goes to — oo as A \ 0. Comparing (|4.20|) and (|4.2ip clearly there exists a / and A > 
sufficiently small such that the corresponding density Uq satisfies (|4.18p . 

Notice that, the greater is M, the greater \a\ has to be chosen for the left inequality in (|4.18p to be 
satisfied. Nevertheless, in case of large mass M, if the two densities are concentrated enough, i.e. if A 
is small enough, the potential in (|4,2ip is strong enough to ensure blow-up of the solution. 

In order to exhibit an example of density satisfying (|4.17p . we follow the same idea as before but 
with the aim to obtain a corresponding energy £ [h,q] with the entropy term dominating the potential 
term. Therefore, we consider the following sequence of densities 



1 N 



<P 



CLi 



A = jvV(2-d) 



(4.22) 



where ip is defined as above and in addition the family of points (ai)i<i<N is choosen so that ai and 
— en belong both to the family. Again, the densities Uq belong to (L 1 n Lz )(M rf ), have mass equal to 
M, L 2 -norm increasing as A \ and the second moment given by 

i N r 

|x| 2 n^(x) dx = M — ^ \ a i\ 2 + A 2 / \z\ 2 (f(z) dz . 



(4.23) 
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Again we assume that A is chosen such that the supports -B(aj,A) of each contribution in (|4.22|) are 
disjoints and we introduce the notation D x (i,j) = dist(-B(aj, A), B(aj, A)). Computing separately each 
contribution of the energy functional, we obtain 

M \og(N\ d ) + I (p(z) log (p(z) dz , 

JR d 

(4.24) 

// "oO^l — -ld^2 n o(y) dxd V 
JjR d xm d \x-y\ d 1 

= 11 Trr^v{ z ')dzdz + 1 1 <p(z)~ t n ,mw y{z')dzdz' 

iV xK / 1 ; |z- Z '| d - 2n N 2 f£ J 7rW \a t - aj + \(z - z')\ d ~ 2 

< I , Viz) z ,, d _M z ') dzdz ' + JH E(^(^')) 2 ~ d » ( 4 - 25 ) 
J JM. d xR d \ z ~ 2 \ JV ~rr. 

due to the choice NX d ~ 2 = 1. 

We claim that there exists a family of points (ai)i<i<Ar such that the last contribution in the r.h.s 
of (|4.25p is uniformly bounded with respect to N. The argumentation goes as follows. First we may 
change the reference norm, i.e. we can replace the euclidean norm | • I2 in W 1 with the supremum 
norm | • 1^, up to some constant. Then, we distribute N points on the regular grid N~ l l d ■ 7L d inside 
the hypercube [— 1/2, l/2] d . Next, we observe that for any index i and any integer k < ^N 1 ^ there 
are at most Ck d ~ l indices j such that | N-^ d k, where C is a constant depending only on 

the dimension. As a matter of fact, after rescaling space by a factor N l / d , those points are regularly 
distributed on a sphere of radius k. Finally, we have the following estimates as N — > +00: A <C N~ l l d 
and consequently D x (i,j) ~ |aj — a,j\i. To conclude, we can estimate the last contribution in the r.h.s 
of (I4T23D as follow 

^£(^)) 2 - d <C^iV £ k^{N-^k) 2 ~ <C±N^\n^ = \c, 

v£j k=l 

where [2~ 1 N 1 ^ d ] denotes the integer part of 2~ 1 N 1 / d . Therefore the interaction potential is bounded 
from below, while the entropy (I4.24|) is decreasing towards —00 as N — > +00. Moreover, it is always 
possible to scale appropriately the location of the family (here inside [—1/2, 1/2] d for the sake of 
reference) to ensure that the right inequality in (|4.17p is satisfied. 

Remark 4.7 (Exotic construction) The above example (|4,22p . which is made of a large superposition 
of approximations of the idenity with disjoint supports, is not likely to be a configuration encountered 
genuinely. We will see later in Section [6] an (over-)simplified version of the parabolic-elliptic (|KSp 
system for which the second blow-up criterion (|4.7p always contains the first one (|4.5p . 



/ 



rio(x) log Uq^x) dx 



1 



N 

v 5 



<p(z) lo§ 



N 



■(f(z) I dz 



and 
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5 Concentration phenomenon for the parabolic-parabolic system 

In this section we shall prove Theorem 1 1.21 i- e - we shall show that the parabolic-parabolic Keller-Segel 
system (e > 0) exhibits a concentration behaviour under some smallness condition on the second 
moment of no, generalization of condition (j4.7|) obtained in the parabolic-elliptic case. However, 
although we are able to obtain blow-up of the cell density in the latter case, we fail to do so as soon as 
e > 0, and we derive a weaker result claiming that the L d l 2 — norm of cell density reaches high levels 
(at least for small e). 

Pretty few results are available for the parabolic-parabolic system (|KSp . We stress out a recent 
contribution by Cieslak and Laurengot [T7] in which they prove the occurrence of such parabolic- 
parabolic blow-up events in one space dimension but with suitable nonlinear diffusion. They follow 
a perturbation strategy, as we do here. However they strongly use the property that the free energy 
functional is bounded from below, which unfortunately does not hold true in our context, at least 
in the usual sense (see Lemma l5.2p . This makes the analysis more difficult and justifies a priori the 
weaker result which is obtained here. 

We consider throughout this section e being positive but a being zero, since our aim here is to 
emphasize the parabolic character of the chemical equation. Then, our key strategy is to consider the 
parabolic-parabolic system as a perturbation of the parabolic-elliptic system. Doing that, the analysis 
of the time evolution of the second moment leads naturally to carefully estimate the time evolution of 
the corrected energy J-[n, c] defined in ()1.6f) . It will be shown that J-[n, c] satisfies a master differential 
inequality (see (|5.12p below), generalization to the parabolic-parabolic Keller-Segel system of the 
differential inequality for J-[n] obtained in the parabolic-elliptic case (see Remark 14. 6p . The master 
inequality together with the second moment growth rate estimate, give us the concentration result. 
Finally, the perturbation strategy could not work without the following two fundamental lemmas. 

Lemma 5.1 (Free energy minimization). Let d > 3. Let n be any nonnegative function in 
(L 1 n L 2d ^ d+2 ">)(M d ), such that / Rd ra(:r) log n(x) dx < oo. Let c := Ed*n. Then, the energy functional 
defined in (jl.4p with a = 0, satisfies 



£[n, c] >£[n,c]= \ 



JR d 



n(x) log n(x) dx 




n(x) c(x) dx 



(5.1) 



for any c such that nc G L l (R d ) and |Vc| G L 2 (R d ). 



Proof. In order to prove (|5.ip formally, it is sufficient to observe that 



[n,c]-£[n,c] = \ I \Vc{x)\ 2 dx - \ [ \Vc(x)\ 2 




n(x)(c(x) — c(x)) dx . 



(5.2) 



Since —Ac = n, a simple integration by part in (|5.2p gives us 
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and the inequality in (|5.1|) is proved, as well as the equality by applying again an integration by 

d 

part. Next, to rigorously justify the previous computation, we observe that Ed G L^~ 2 (M ) and 

\VE d \ G Lw (R ). Therefore, we have on the one hand c £ L d - 2 (K ), and on the other hand 
|Vc| G L 2 (R rf ) by the weak Young's inequality. Consequently, the integration by part in (|5.2j) can be 
performed (see for instance [39]). □ 

Lemma 5.2 (Corrected energy lower bound). Let n be any nonnegative function in (L l C\L?- )(M. d ) with 
finite second momentum. Let c be such that nc G L 1 (M. d ) and |Vc| G L 2 (K rf ). Wii/i M = ntt^) 
and i/ie constant B(d,M) defined in (|4.12p . the following lower bound for the corrected energy ^[n^c] 
holds true 

d(d - 2)MT[n, c] + B(d, M) + tJz^M C HLS (d, d - 2) \\n\\ L a > 2dM . (5.3) 



Proof. Let c = E d *n. Using the free energy minimization Lemma [5. II and the expression for B(d, M), 
we have 



d{d-2)MJ r [n,c] +B(d,M) > 2(d - 2 

+ 2(d-2)£[n,c 
= 2(d-2) 



dM, r „„ „, dAf, (dM 

— log J + — - M log M - — log I — 



+ 2dM 



f , r ^ ,r dM , fdM\ 

J n log n + —(log / + 1) - M log M — log f — 1 

+2dM - (d- 2) / nc. (5.4) 



Since the bracket in the r.h.s. of (|5.4p is nonnegative by the entropy lower bound (|2.9p . inequality (|5.3p 
follows from (|5.4f) when applying to the quasi-stationary potential f^ d nc, the potential confinement 
inequality (pOOD . 

□ 

Before proceeding in proving Theorem II .2^ let us observe that the constant C(d) in (jl.9p is the 
same constant as in the master inequality (|5.12p . It derives from explicit computations in the proof 
of Theorem 11.21 Therefore, it can be quantified. Moreover, it is important to underline that the 
initial concentration condition (jl.8p from one hand converges to the blow-up condition (|4.7p as e — > 
and from the other hand it implies that e is small with respect to ||no||^d m<JN . Indeed, applying the 



logarithmic function to both side of (jl.8p . after rearranging the term, (jl.8p reads as 

d{d- 2)MT[n ,c ] +B(d,M) + d(d - 2)Me 7 < 0. (5.5) 
Then, the corrected energy lower bound (|5.3p gives 

d(d-2) £ T + 2d< \$ d - 1 \- 1 C H Ls(d,d-2)\\no\\ d (5.6) 
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Obviously, inequality (|5.6|) does not imply that Theorem 11.21 holds true only for small s > 0. Never- 
theless, inequality (jl.9p gives a relevant concentration result only if e is small and ||no|| d is upper 



'£2( 

bounded by the constant given in the r.h.s. of (|1.9|) . One can try to obtain a concentration result 
for n(t) when e is large, for exemple taking 7 > 1 or a more general function of e instead of e 7 into 
(jl.8p . However, the only result obtained is that the L d / 2 -norm of n(t) stays large for a while when 
the L d l 2 -norm of no is large (see also Remark 15. 3p . Finally, the maximal time of existence T max in 
Theorem 11.21 is not requested to be finite and we cannot exclude that the solution (n,c) is global in 
time. 

Proof of Theorem M.SX We shall proceed in several steps to conclude with a contradiction argument. 
The delicate and original points of the proof consist in the perturbation argument of the first step and 
in controlling the growth of the second moment I{t), (no larger than ilogt), in the third step. 

First step : the chemical quasi- stationary approximation. We express the chemical c as 

c = E d *{n-ed t c). (5.7) 
Indeed, applying the Fourier transform to the parabolic equation for c, we obtain 

c(6 = ^(n-^c)(0. 

Then, it suffices to apply the inverse Fourier transform and Theorem 5.9 or Corollary 5.10 in [39], to 
obtain the identity (|5.7|) . Consequently, the gradient of c can be written as follows: 

Vc(M) = -|Sg=iT / , X ~ y ,A n(y,t)-ed t c(y,t))dy. (5.8) 

l§ I jR d \ x - y\ 



Second step : growth o/J-"[n,c]. Introducing the representation (15. 8p into the evolution equation (14.2 
for the second moment I of n, we get after symmetrization 

^J(t) = 2dM--? zr [[ n(x,t) X '^ X ~ f n(y,t)dxdy 
at \S d L \ JJ R d xR d \x - y\ d 

\ I \ * n(x ' l) X \} x Jd 9 * c ^' *) dxd y 

P I J JR d xR d \ x — y\ 



2dM - (d - 2) / n{t)c{t) 



+T^Tu II n(M) , ,? d t c(y,t)dxdy. (5.9) 
P I J jR d xR d \x-yr 



where c = Ed * n. In order to control the second integral term in the r.h.s. of (15 .9p . we apply the 
general Hardy-Littlewod-Sobolev inequality |39| . and we use a suitable Holder type inequality to get 
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for every 5 > 0, 

// \x\n(x,t)- \d t c(y,t)\dxdy < C(d)\\xn(t,x)\\ ^ \\d t c{t)\\ L2{Rd) 

JjRd xR d \X-y\ d 1 Ld+2( M d) v > 



< C(d)l2(t)\\n 



\d t c 



L 2 



< C(d) 



s IK*)IU^ + s~ 1 i(t)\\d tC (t)\\ 



L- 



,(5.10) 



where the constant C(d) above may change from line to line, the last being the one appearing in (jl.9p . 
Plugging estimation (|5.10p into (|5.9p . we obtain our first main estimate (the equivalent of (|4.3p for 

e > 0) 

j t I(t) < 2dM - (d - 2) n(t)c(t) + e5C(d) \\n(t)\\ ^ + e5~ 1 C(d) I (t)\\d t c(t)\\ 2 L2 {Rd) . (5.11) 



Next, we use the free energy minimization Lemma 15. II and the free energy dissipation equation ([1.5 
into (]5.1ip to have 



d 
dt 



I(t) < 2dM+2(d-2) 



£[n,c](t) - / n(t) log n(t) 



+s5C(d)\\n(t)\\ d{Rd -5- 1 C(d)I(t) j t £[n,c](t) 



Choosing 5 := — C(d), the entropy lower bound (|2.9|) gives us 



—I(t) < 2dM + 2(d- 2)£[n,c]{t) + 2(d - 2) 



dM, r dM , r 
— log I + — - M log M — log 



dM , /dM\ 



i.e., after rearranging the terms, we get the following differential inequality for the corrected energy 



+ -dMC 2 {d)\\n{t)\\ „ „ --Lm)4f[n,cim 



J(t) — T[n,c){t) < d(d - 2)M T[n,c](t) + B(d,M) + ^-dMC 2 {d) \\n\z 



(5.12) 



Third step : moment's growth. Here we compute the evolution of I{t) in an alternative way, involving 
the Fisher's information in the the energy dissipation equation (jl.5p . As a consequence, we shall 
prove that the second moment does not increase asymptotically faster than 2dMtlogt as long as 
||^^ WN stays bounded from above. Indeed, 



dt 



lit) 



< 2 



n(x, t) x ■ V(log n(x, t) — c(x, t)) dx 

n(x, t)\ V(log n(x, t) — c(x, t))\' z dx 



\x\ 2 n 



\ 1/2 / r \ 1/2 

(x, t) dx I 



Therefore, 



il 1/2 (t) <( [ n(x,t)\V (log n(x,t 
dt \Jud 



) -c(x,t))\ 2 dx 



1/2 
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Integrating the above inequality over (0, t) we get 

1/2 



and 



I 1/2 {t) < I 1/2 (0) + (^j n{x,s)\X7 {log n{x,s) - c(x, s))\ 2 dx^J ds . 
I(t) < 21(0) + 2t [ [ n(x,s)\V (log n(x,s) - c(x , s))\ 2 dxds . 

JO JR d 



Using the dissipation of the energy, we derive the following pointwise estimate for I(t) with respect 
to the corrected energy J-"[n, c] (t) 

I(t) < 21(0) + 2t(£[n Q ,co] -£[n,c](t)) = 2/(0) + 2t£[n Q , c ] - dMt T[n, c] (t) + dMt logl(t). 

Finally, being dMt log I < \ l + dMt log(2dMt), we obtain 

I(t) < 4/(0) + 4t£[n ,co] - 2dMt T[n,c](t) + 2dMt \og(2dMt) , 

i.e. the claimed behaviour for I(t), thanks to the corrected energy lower bound ()5.3fl . 



I(t) < 4I(0)+4t£[n ,c ]+2t 



(d - 2)~ l B(d, M)+fi d M C HLS (d, d - 2) \\n\z) \\ T e. 



+2dMt log(2dM t) . 
(5.13) 

Fourth step : conclusion. We now conclude showing by a contradiction argument that the concen- 
tration result (ll.9p holds true. Indeed, comparing the master equation (|5.12p and the concentration 
condition (11.81) written as (|5.5p . we claim that we can not have uniformly in time 

£ -dMC 2 {d)\\n{t)\\ Li <d{d-2)Mei . (5.14) 



In order to proceed, we have to distinguish between two cases according to the fact that the initial 
density no satisfies or not (|5. 14j) . If no does not satisfy (|5.14p . the concentration result (jl.9p is 
obvious. If not we deduce from the master equation (|5.12|) and from (|5.5|) that the corrected energy 
is initially decreasing. Next, if (|5.14p holds true uniformly in time, there exists 5 > such that 
/(t)^^ 7 ^, c](t) < —5, uniformly in time, i.e. F[n, c] remains decreasing for t > 0. Moreover, plugging 
the upper bound for ll^Wll^ce,,^ given by (|5.14p into (|5.13p . we get the estimate 



d S 

^[n, c](t) < - 4/(0) + C ( £) d,M,£[n ,c }) t + 2dMtlogt ' (5 ' 15) 

Since the right-hand-side of (|5.15p is not integrable at infinity, but J 7 [n,c](t) is bounded from below 
as soon as ||n(t)||^d m<n is bounded from above. We have obtained a contradiction which completes 



the proof of the theorem. □ 
Remark 5.3 It is worth noticing that when no satisfies (|5.14p . from (|5.6p e results necessarily small, 

£< d\§*-i\c*(d) • (5 - 16) 
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Moreover, the above upper bound for e is independent from the function of e choosen in the initial 
concentration condition (jl.8p . On the other hand, if uq does not satisfy (|5.14|) . we have 



Tin d 
1 U|I L^ 



> max \d{d - 2)\S d - 1 \C„ 1 LS (d, d - 2) ; e^ 1 1 



and the case "e large" is included here. 

Remark 5.4 (Concentration result of order e" 1 ) It is possible to obtain a concentration result of 
order e _1 under the "critical" hypothesis 

d{d-2)MF[n Q ,c G ] + B{d,M) + £ -dMC 2 {d)\\n4 L d m <0, (5.17) 

that writes equivalently as 

\x\ 2 n (x)dx < K 2 (d) M 1+ d exp ^ _ ^^[ n o, c ]^ exp ^" ^(^ll^oll^j 

Then, from the master equation (j5. 12[) we obtain that the corrected energy J-[n, c] is initially decreas- 
ing. Proceeding as before by contradiction, we deduce the existence of a time T > such that 

d(d - 2)MT[n , c ] + B(d, M) + | dMC 2 {d) \HT)\\^ > 



However, the critical concentration initial condition (|5.17p implies that e is necessarily small, i.e. e 
must satisfies (|5.16p . 

Remark 5.5 Theorem 1 1 . 2 1 does not give any clue about the solution's behaviour after the first "con- 
centration time" (blow-up, persistence, dispersion ?). This is due to the fact that the master equa- 
tion (|5.12p does not bring any information about the evolution of the second moment I(t) and of 
||n(t)|| d after that time. 



Remark 5.6 It would be possible to extend the first blow-up criterion given in Proposition 14.11 to 
the fully parabolic system as well. However the method is less natural and the final result is weaker. 
Thus we will not develop the arguments in this paper. 



6 A discrete model mimicking the parabolic-elliptic system 

In this section we aim to give a geometric intuition of what could be the dynamics of the parabolic- 
elliptic Keller-Segel system in space dimension d > 3. For this purpose and for numerical convenience, 
first we shall replace the classical Keller-Segel system (|KSp with e = a = 0, by a one-space dimensional 
variant having an analogous behaviour to (]KSp . Energy and especially homogeneity considerations 
have guided our choice for the substitutive system as we shall see below. Next, we shall rephrase the 
one dimensional Keller-Segel system using the pseudo-inverse distribution function of the cell density n 
(see [3 [26]). A numerical scheme following an idea of [7] is discussed and a discrete dynamical system 
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is then proposed, with the aim of providing a visualization tool for the dynamics of the Keller-Segel 
system. 

Let us recall the action of the mass-preserving dilation f(x) — > f x (x) = \~ d f(\~ 1 x) on the free 
energy (jl.7p with a = 0, yet used in the subsection 14.21 

£[n x ] = — dM log A + j n(x,t)logn(x,t) dx — ^— — ff n(x,t)Ed(x — y)n(y,t) dxdy . (6.1) 

We clearly see the different homogeneities of the two contributions composing the energy, namely the 
entropy and the potential. The first one has say an "almost zero" homogeneity (up to a logarithmic 
correction), while the second has the same degree of homogeneity as the kernel Ed- These different 
homogeinities make the behaviour of the (|KSp system more intricated in high dimension than in 
2 dimension (see Remark 16. 2p . In order to reproduce the same high-dimensional behaviour in the 
simpler frame of the one space dimension, we consider from now on the following one parameter 
family of systems describing self-attracting diffusive particles 

I d t n = d xx n - d x (nd x c) , . 
c = K-f * n , 

where the interaction kernel is given by K~(x) := 7 -1 |^| _7 , with 7 G (0, 1). As for the (|KSp system, 
the new system (|6.2p is equipped with the decreasing free energy 



£j[n](t) = / n(x, t) log n(x, t) dx — - / n(x,t) (K y * n)(x,t) dx , (6.3) 
Jr 2 7r 

satisfying 

d f 

— £j[n](t) = — / n(x,t) \d x (log n(x,t) — (K 7 * n)(x,t))\ 2 dx . 
dt 7r 

The action of the one dimensional mass-preserving dilation f(x) — > f x (x) = A" 1 f(X~ 1 x) on the new 
energy (|6.3p is then given by 

£[n x ](t) = —M log A + f n(x,t)logn(x,t)dx-^— f n(x,t)(K^ * n)(x,t) dx , (6.4) 

Jm. 2 7r 

so reproducing (|6.ip . 



Remark 6.1 The case 7 = —1 is not included here due to integrability issues. It is too singular in 
dimension 1, although it would naturally correspond to the Poisson kernel in dimension d = 3 as far 
as homogeneity is concerned. The generalized model (j6.2[) will be the subject of another paper [12j 
(not restricted to dimension 1 only). 

Remark 6.2 (The two dimensional case) Let us recall that in the two dimensional case the interaction 
kernel (jl.2p for the parabolic-elliptic (JKS]) system with a = is of logarithmic type, i.e. E2(x) = 
— ^ log \x\, and the energy (jl.7p reads as 

£[n](t) = / n(x,t)logn(x,t) dx + — I n(x, t) log \x — y\ n(y,t) dx dy . (6-5) 
7r2 47r77 R 2 xK 2 
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Therefore, under the two dimensional mass-preserving dilation, the two energy contributions have the 
same "almost zero" homogeneity (up to a logarithmic correction) and 

£ [n x ] = 2M (¥- - lj log A + £[n] . 

The above property, in the gradient flow interpratation, directly implies blow-up of any density n 
having super-critical mass M > 8ir (see [71 112]). 

Remark 6.3 (The limiting case 7 = 0) It is possible to mimic the above two dimensional case with 
the help of the one dimensional system (|6.2p . if one replace the kernel K 7 with K 1 (x) = 7 _1 (|x|~ 7 — 1) 
and let 7 go to 0. Doing that, the limit interaction kernel is given by Kq(x) = — log | x | and the 
corresponding free energy satisfies 

£0 [n X ] = M - lj log A + So [n] , (6.6) 

under the one dimensional mass-preserving dilation, exactly as for the classical two dimensional Keller- 
Segel system (|KSp , As a consequence, system (|6.2p . with this particular choice of kernel, yields a critical 
mass phenomenon, M = 2 being the critical threshold, as we shall see below (see Remark I6.6|) . This 
system has been considered previously in [6] I14j. 

Remark 6.4 (Corrected energy homogeneity) It is worth noticing that, still in the parabolic-elliptic 
case with a = 0, the corrected energy J-[n] given in (|1.6p . has the following striking homogeneity 
structure: the action of the mass-preserving dilation on the term log ( L d \x\ 2 n(x,t) dxj cancels with 
the logarithmic entropy's contribution, so that 

F[n x ]= I n(x,t) log n(x,t)dx - * [[ n(x,t)E d (x - y)n(y,t) dxdy . 
JR d dM J J R d xR d 

6.1 Gradient flow interpretation of the one dimensional system 

Let us define the pseudo-inverse distribution function of any nonnegative L 1 (M) function / having 
mass M, as the function X(m) defined on the interval (0, M) as follows 



f(y)dy>m\, me(0,M) 



X(m) := inf lx € I 

Then, for the pseudo-inverse distribution function X = X(m, t) of the density n, the energy functional 
(16.31) rewrites as 



i-M 1 rM pM 

£ [n](t) = g 7 [X](t) = - log (d m X(m,t)) dm / / \X(m,t) - X(m f ,t)\^ dmdm' , (6.7) 

Jo 2 7 Jo Jo 

while system (|6,2p rewrites as follows [26] : 

-d t X(m, t) = d m ( v ) ) + [ M sign (X(m, t) - X(m', t)) \X(m, t) -X(rri, t)!" 7 " 1 dm' , (6.8) 
\d m X(m,t)J Jo 
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Furthermore, the action of the one dimensional mass-preserving dilation on the pseudo-inverse dis- 
tribution function simply consists in multiplying by A, i.e. X x (m) = XX{m). Therefore, the same 
different homogeneities of the two contributions in the new energy functional (|6. T[) clearly appear 
again 

pM \—f pM pM 

7 [X A ](i) = -MlogA- / log (d m X(m, t)) dm- — — / / \X(m, t) - X(rri, t)\~~< dm dm' . 

Jo 27 J J 

Proposition 6.5 (Gradient flow interpretation). The integro- differential equation (|6.8p is the gradient 
flow of the energy functional Gj[X] for the Hilbertian structure on L 2 (0, M): 

d t X = -V L 2g 7 [X] . (6.9) 

In addition, the L 2 (0,M) norm of the pseudo-inverse distribution function X of n is the second 
momentum of n: 

ll^(-,*)IU 2 ((W) = / \x\ 2 n(x,t)dx . 



This interpretation is due to Jordan, Kinderlehrer and Otto |33l 144] . We refer to the books [13 [T] 
for a comprehensive presentation of this interpretation in any dimension of space. Here, for sake of 
completeness, we give a sketch of the proof of identity (j6.9f) . 



Proof. To prove (|6.9p . we compute formally the variation of the functional directly from ([67 

to obtain 

r-JW 



PM 

g i [X + H] = - log (d m X(m) + d m H(m)) 
Jo 



dm 

-1 pM pM 

- / \X(m) - X(m') + Him) - H(m')\-~i dmdm' 
27 Jo Jo 



M f M 



+ U I |x(m) - xim ' r WfM dmdm ' + 0{H) 



/■M pM H(m) 

+ L L \ x{m) - x{m ' ) \~ \( m) -k m ') imdm ' +0(H) - 



□ 



Remark 6.6 (The critical mass phenomenon in the limiting case 7 = 0) The energy functional C?o[-X]j 
given by (|6.7p when replacing the kernel with Kq = — log \ x\, satisfies 

g [xx] = m(^-- i\ log a + g [x] , 
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as one can deduce for exemple from (|6.6p . Differentiating the above relation with respect to A and 
evaluating the result at A = 1 yields 



(x,vg 7 [x]) i2(0M) =M(^-i) 



On the other hand, (X, V^ 7 [X])£,2( 0) m) is nothing but the time derivative of \\X(-, i) || z, 2 (o,JW) ( U P 
to a change of sign), as shown in Proposition 16.51 Therefore, we recover that the positive function 
||-X"(">£)IIl 2 (o,a/) has a linear time decay when mass is super-critical, i.e. M > 2. 

6.2 The steepest descent scheme and the simplified dynamical system 

In [33] the authors introduce a steepest descent numerical scheme for the linear Fokker-Planck equa- 
tion. This coincides with the Euler time-implicit scheme when the equation is reformulated using the 
pseudo-inverse distribution function. Those schemes possess the advantage of capturing important 
energy features (energy is decreasing and uniform in time estimates can be derived). 

This strategy has been adapted to systems of Keller-Segel type in [7J, and particularly to the one- 
dimensional variant (|6.2p with the logarithmic interaction kernel Kq. It consists in discretizing the 
energy on a regular mesh in the mass space (0, M), and then performing a gradient flow for the finite- 
dimensional functional obtained. As the numerical scheme is carefully built to preserve the gradient 
flow geometry of the original system, this procedure may give a clear intuition of the problem. In the 
sequel, we shall construct the same type of numerical scheme but starting from the energy functional 
(16771) with 7 > 0. 

For sake of simplicity, we opt for a very rough discretization of the space (0, M), namely we consider 
the following three points regular mesh 

M M 3 

mi = — , m 2 = — , m 3 = -M . 

We next discretize the energy using finite differences. This writes, up to a constant factor h, as 

G 7 [Xi, X 2 , X 3 ] = - log(X 2 - X{) - log(X 3 -X 2 )-- {{X-, - X 2 )~~i + {X 3 - + (X 2 - X^) , 

(6.10) 

where h = M/4 is the space step, while the gradient flow of (|6.10p is given by the discrete dynamical 
system 

( = ~ x 2 (t)-x 1 (t) + h ((XsW " Xl (t)) " 7 " 1 + {X2[t) " Xl ^))~ 7 " 1 ) 

X ^ = x 2{ t)-x l{ t) - m)-x 2{ t) + h ^ ~ x > w)" 7 " 1 - ^ w - x ^~^ 



Xs(t) = Xi{t) l _ x<t) ~ h ((^3(t) - X 2 {t))-^ 1 + (X 3 (t) - XiCt))^" 1 ) 



(6.11) 
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Let us observe that everything can be rewritten in terms of u(t) = X 2 (t) — X\(t) and v(t) = 
X$(t) — X2(t) because the center of mass X\(t) + X2(t) + X^(t) is conserved and assumed to be zero 
here, without loss of generality. Furthermore, we can reconstruct X±(t), X 2 (t), X$(t) from u(t) and 
v(t) by the following relations 

Xl (t) = ~(2u(t) + v(t)) , X 2 {t) = i (u(t) -v(t)) , X 3 (t) = ^(u(t) + 2v(t)) . 



(6.12) 



The dynamics of the new variables (u(t),v(t)) are given by 

= W)~W) +h ^ (t) " 7_1 " 2 ^)~ 7_1 - («(*) + ^))~ 7_1 ) 

= ^ - + h (u(tr^ - 2v(tr^ - (u(t) + v(t))-^) 

However, system fj6. 12f) is not a gradient flow in the updated variables (u, v). 

As an example we plot in Fig. [I] the case 7 = 1/2 and M = 1.6 (it would be similar for other choices 
of M > and < 7 < 1). One clearly observes the separation between two opposite alternatives: 
"blow-up" when u — > or v — > and "dispersion" when both u — > +00 and v — > +00. 

6.3 Investigation of the landscape's geometry: blow-up vs. global existence 

In this section we shall derive two blow-up criteria and one global existence result for the dynamical 
system (16. IIP . For that, it will be useful to rewrite the squared euclidean norm of X = (Xi,X 2 ,X 3 ) 
in term of the updated variables (u,v), 

\X\ 2 = - (u 2 + v 2 + uv) = - (u 2 + v 2 + (u + v) 2 ) , (6.13) 
3 3 

as well as the functional G~[X] 

G y [X] = -logu-logu- ^ (u-^' + v-^ + (u + v)-t) , X=-£- (6-14) 

Proposition 6.7 (Blow-up criteria). Assume that the initial point Xq satisfies one of the two following 
criteria: 

|X ° |2< (f) 2/7 ' (6 - 15) 



\X \ 2 < 2exp I -G 7 [X ] - - , (6.16) 
then either u(t) or v(t) vanishes in finite time. 

Proof. We follow carefully the evolution of the euclidean norm |X(t)|. The functional G 7 can be 
decomposed according to homogeneities: 



G^[X]=V[X]-W y [X], < 



XJ[X] = - log(X 2 - Xi) - log(X 3 - X 2 ) , 

y (6-17) 
W 7 [X] = ^ ((X 3 - X 2 )-> + (X 3 - + (X 2 - Xi)" 7 ) . 
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Figure 1: Behaviour of the dynamical system f|6.12[) in the phase plane (u, v). (Top) The first f (|6.15j) ; left) resp. 
the second ( (|6.16[) ; right) blow-up criterion is figured as a bold line. (Middle Left) The line having equation 
W 7 [X] = f (where the functional takes maximal values along radial rays) is plotted in bold. It is postulated to 
be an admissible criterion for blow-up but this has yet to be proved (see (|6.27|l and the following discussion). 
(Middle Right) The global existence criterion is figured as a bold line. (Bottom Left) The unstable manifold 
(bold line) starting from the energy maximal point clearly separates the two basins of attraction. The derivation 
of its equation would provide a single criterion to distinguish between global existence and blow-up. (Bottom 
Right) Finally all the previous lines are plotted on the same figure for the sake of comparison. 
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We derive the following Euler formula under mixed homogeneities assumptions: 

<X,VG 7 [X]) = Ag 7 [AX] u=1 = - 2 + 7 W 7 [X] . (6.18) 
Therefore, the evolution of the euclidean norm under the gradient flow is driven by: 

\j t \X{t)\ 2 = - (X(t), VG 7 [X(i)]) = 2 - 7 W 7 [X(t)] . (6.19) 
We use the following Jensen inequality based on the convexity of the function (•) ~ 7 / 2 : 

Q {u 2 +v 2 + {u + v) 2 )^j V < i( u -T + t ,-7 + ( u + u )-7) , (6.20) 



From (|6.19p we obtain: 



11 "3x 7L J 



±±\X(t)\ 2 <2-3 X \X(t)r. (6.21) 



Therefore, the norm necessarily vanishes in finite time if the initial data satisfies the following criterion: 
\X \ < (3x/2)V7. 

Starting from equation (|6.19|) we can plug the energy in the computation: 

\j t \X(t)\ 2 < 2 + 7 G 7 [X ] - jU[X(t)] . (6.22) 
We use an alternative Jensen inequality, based on the concavity of log(-): 

log Q (u 2 + v 2 + uv)^j > log 2 + i(21ogu + 21ogu + log(m;)) , (6.23) 
log|X| 2 > log2-U[X]. 

We obtain consequently: 

\j t \X{t)\ 2 < 2 - 7 log 2 + 7 G 7 [X ] + 7 log \X(t)\ 2 . (6.24) 

Therefore, the norm necessarily vanishes in finite time if \Xq\ 2 < 2exp(— G 7 [Xo] — 2/ 7 ). □ 

Remark 6.8 (Comparison between the two criteria) As opposed to Section the two criteria (16. 151) 
and (|6.16p are not complementary. Indeed the former implies the later (see Fig. [TJ. We shall prove 
this in the sequel. However there is a first point when looking for equality cases in resp. (|6.20p and 
(|6.23p . In fact, the former admits not equality case, whereas the later is an equality when u = v. 

To prove that the first criterion (|6.15p enhances the second (|6,16p . we shall prove more generally 
that for all Xq, 

G 7 [X ] + H < log 2 - ^log . (6.25) 
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_7' 1 1 1 1 ' 

0.2 0.4 0.6 0.8 1 

u = Simplex parametrization u = 1 

Figure 2: Gauge function for (resp.) the three-points and four-points discretization of the energy functional. 
Recall that it is a function of homogeneity zero, which is defined by (|6.27j) . where TJ[X] = — logu — logi; (three 
points) or TJ[X] = — logu — logv — \ogw (four points), and W 7 [X] = X7 1 ("~ 7 + v 1 + (u + v)~~<) (three 
points), and W 7 [X] = X7 _1 ( u ~ 7 + w ~ 7 + w ~ 7 + ( u + + ( v + w)~~' + (u + v + w)~~ < ) (four points). The 
maximal value is plotted on the simplex (resp. u + v = 1, u + v + w = 1) because it is a functional of homogeneity 
zero. 



The maximum of G 7 is achieved on the diagonal {u = v}. A simple computation yields for the 
extremal point: 

-7 2 

* x(2 + 2-7)' 
Plugging this into (|6.25[) we are reduced to prove: 

The function 7 1— > log (2 + 2~ 7 ) being convex, we have for 7 > 0: 

log(2 + 2- 7 ) -log3> -|log2> -|log2. 

To conclude, the second criterion involving energy is clearly better in this case. This strongly uses the 
fact that the free energy is bounded from above (together with a precise evaluation of the maximal 
value), whereas this is no longer true in infinite-dimension. These computations are possible due to 
the simplicity of the toy model. 

Remark 6.9 (Derivation of a single criteria: open issue) The separation between the basins of at- 
traction of the two axes (blow-up) and of 00 (global existence) is given by the unstable manifold of 
the energy critical point drawn on Fig. [TJ We have failed in deriving an equation for this manifold. 

Remark 6.10 (Towards a better description of the landscape's geometry) We introduce below another 
correction of the free energy. Notice that although it comes out very naturally, it is still unclear 
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how to benefit from it appropriately. We aim to gain some insights concerning the geometry of 
the functional landscape, by means of homogeneity arguments. We decompose the free energy as 
previously: G 7 [X] = U[X] — W 7 [X] where U is almost zero- homogeneous (up to a logarithmic 
correction) and W 7 [X] is (— 7)-homogeneous. We define the maximal value along rays: 

JLy[X] = maxG 7 [AX] = TJ[X] - - (log h W ^ \ + i 

A>o 7 \ V 2 / 

This function is zero-homogeneous by definition. Moreover, a simple computation shows that the 
set of extremal points: {X*X \ X* = argmax^ G 7 [AX]} admits the following equation: W 7 pT] = 2/7 
(cf. (|6.18p ). Observe that the maximal value coincides with U[X] — 2/7 precisely on this curve 
(hypersurface in higher dimension). 

A reasonable claim would be that the gradient flow can not overpass this "maximal line", so that 
the later splits the phase space into two regions, and blow-up is guaranteed whenever we start initially 
on one side of this line. It turns out that this "maximal line" does not coincide with the unstable 
manifold which is the clear separation between global existence and blow-up (see Fig. [1]). 

It is also worth noticing that the set {W 7 [X] > 2/7} coincides exactly with the area of the phase 
space where the second momentum is initially decreasing (|6.19p . The connection obviously derives 
from the definition of H 7 [X] . 

To conclude we give a global existence result which is somehow analogous to Theorem 11.11 This 
criterion is optimal in the sense that the asymptotes coincide with the unstable manifold (not shown). 

Proposition 6.11 (Global existence). If the initial data satisfies: 

X ( V + V) < 1 , (6-28) 

then the solution is global. 

Proof. To prove global existence in the line of Section [3] we shall estimate some suitable quantity pre- 
venting u(t) and v(t) from vanishing. Namely, we investigate the evolution of Lit) = 7 -1 (u(i)~ 7 + f(i)~ 7 
We have: 

j t L(t) = - [2u{t)-^- 2 + 2v{t)~^ 2 - u(t)- 1 u(t)-T- 1 - u{t)- x v(t)-*i- x ] 

+ X [2u(t)~ 2 ^ 2 + 2v(t)~ 2 ^ 2 - 2v{t)-"'- l u{t)-"i- 1 + (u(t) + vit))-"/- 1 (u(t)" 7 " 1 + v(t)-^ 

We claim that the following inequality (16.29P holds true. This enables to conclude that the quantity 
L(t) is nonincreasing in time if (|6.28p is verified. 

2u~ 2 ~ / - 2 + 2v~ 2 ~ / - 2 - 2v^~ 1 u~^ 1 + (u + w)" 7 " 1 (n" 7 " 1 + v'^ 1 ) 

< [-u" 7 + v~~<] [2U" 7 " 2 + 2t>- 7 " 2 - v^u^- 1 - u^v^- 1 ] . (6.29) 



(6.27) 
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This is obvious when homogeneities coincides (7 = 0). Due to homogeneity, (|6.29[) reduces to the 
following inequality for U = u/v: 

2t/ ._a 7 -2 +2 _2 t r7-i + (£r + i)-7-i (U-1- 1 + 1) < [C/~ 7 + l] [2C/~ 7 ~ 2 + 2 - U^' 1 - U' 1 ] . (6.30) 

We restrict to the case U < 1 without loss of generality. We bound the delicate term (U + 1) _7_1 by 
1. We end up with: 

2C /-27-2 + 3 _ j 

on the one hand, and 

2 {7-2y-2 + 2 + 2C/" 7 " 2 + 2J7~ 7 - U' 2 ^- 1 - U~ l - 2[/~ 7 ~ 1 , 

on the other hand. Thus, we are reduced to prove that: 

1 + C/" 27 " 1 + U~ l + U^ 1 < 2C/- 7 - 2 + 2C/" 7 , (6.31) 
1 + f/" 7 - 1 < U-" 1 + U~^ 2 + I/' 1 (C/- 7 " 1 - 1) (1 - U 1 -" 1 ) (6.32) 

the later being clear term by term because 7 € (0, 1) and U < 1. 

□ 

7 Concluding remarks and open problems 

In this paper we have investigated existence and blow-up issues concerning the Keller-Segel system 
(|KS|) in high space dimension (namely d > 3). Concerning the blow-up, we have analyzed two sufficient 
criteria for the parabolic-elliptic Keller-Segel system, and we have discussed complementarity between 
them. We have proposed a possible extension of one of these criteria to the parabolic-parabolic 
case. Although we cannot still exhibit blow-up in that case, we are able to prove that the "critical" 

L d/2 

—norm of the density n reaches high level after some finite time. We have also proposed a discrete 
model for visualizing the dynamics of the parabolic-elliptic Keller-Segel model in finite dimension. 
All this analysis highlights the complexity of the high dimensional Keller-Segel system. Although in 
dimension 2 the existence and blow-up issues, in a simplified scheme, reduce to the mass being smaller 
or greater than 8ir, here it appears that the separation between global existence and blow-up is not 
so easy to describe, even for the discrete model. 

We aim to conclude the paper listing below some problems which are still open up to our best 
knowledge. 

(a) There is a gap between the global existence criterion ()3.30p and any blow-up criteria (|4.5p and (|4.7p . 
Indeed, both of them imply the inequality 

II 1 - 1I z2My)dxdy>2dM\S d ~ 1 \. (7.1) 

J jR d xm d \ x — y\ 
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which implies 

2d|S rf - 1 | 2d 8 

l|no|l ^(M d ) > C HL s{d,d-2) ~ (d - 2) Cf(d) > dC 2 s {d) ' 

Let us observe that (|7. 1 j) reduces to M > 8ir in dimension 2. 

(b) The two blow-up criteria share the following feature: they imply that ^/(i)| t=0 < 0, that is 
nothing else than (|7.ip . Would it be possible to prove that the initial decreasing behaviour of the 
second moment is sufficient for blow-up ? This seems to be the case on the finite dimensional 
model. However this does not give an exhaustive description of the set of initial data for which 
blow-up occurs. 

(c) It has been shown that none of the two blow-up criteria implies the other. However in the discrete 
model, the criterion involving energy is clearly better. Is it possible to show that the energy 
criterion (|4.7|) . which is more convenient to deal with, is also better in some sense for the (|KSp 
system ? 

(d) Are we able to derive the equation for the separation line (unstable manifold) in Figure [Tp Observe 
that such an unstable manifold suggests the existence of a stationary state, at least in the discrete 
model for which the free energy is bounded from above. It is a challenging issue to understand 
better the geometry of the energy landscape for the parabolic-elliptic Keller-Segel system. 
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